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In this paper we address several new developments in the theory of dc Josephson effect in 
superconducting weak links. We analyze an interplay between quantum interference effects and 
Andreev reflection in SNS nanojunctions with insulating barriers and demonstrate that such effects 
may qualitatively modify the Josephson current in such structures. We also investigate an impact 
of the parity effect on persistent currents in superconducting nanorings interrupted by a quantum 
point contact (QPC). In the limit of zero temperature and for the odd number of electrons in the 
ring we predict complete suppression of the supercurrent across QPC with one conducting mode. In 
nanorings with SNS junctions a 7r-state can occur for the odd number of electrons. Changing this 
number from even to odd yields spontaneous supercurrent in the ground state of such rings without 
any externally applied magnetic flux. 
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In 1926 Albert Einstein posed a remarkable question 
[1] : "Of particular interest is the question whether a link 
between two superconductors also turns superconduct- 
ing" . The answer to this question was provided by Brian 
Josephson in 1962 [2]. It was predicted by Josephson 
that dissipativeless flow of Cooper pairs between two dif- 
ferent superconductors separated by an insulating barrier 
is possible provided this supercurrent Is does not exceed 
some critical value Iq- Furthermore, the dependence of 
this current on the phases of macroscopically coherent 
wave functions of Cooper pairs was established in a very 
simple form [2] 



Is = Ic sin ip, 



(1) 



where p is the difference between the phases of the BCS 
order parameters of two superconductors. Eq. (1) rep- 
resents the dc Josephson effect. The answer to the Ein- 
stein's question [1] turned out to be positive. 

What if the total current / flowing through the barrier 
is larger than Ip? In this case a part of the net cur- 
rent across the barrier is transferred by normal electrons 
(quasiparticles) and the rest of it is carried by Cooper 
pairs. While the second contribution, Is, remains dis- 
sipativeless and is again described by Eq. (1), the first 
- quasiparticle - contribution to the current is dissipa- 
tive and, hence, causes a non-zero voltage drop V across 
the insulating barrier. In the presence of this voltage the 
coherent phase difference tp acquires a time dependence 
described by another famous Josephson relation 



dip 
~dt 



2eV 



(2) 



Combining Eqs. (1) and (2) one immediately arrives at 
the conclusion that for any non-zero V the supercurrent 
Is changes in time. In the case of time-independent volt- 
ages one has ip = 2eVt/h and, hence, the Josephson 
current (1) will oscillate in time with the fundamental 
frequency proportional to the voltage V. Eq. (2) and re- 
lated to it oscillations of the supercurrent represent the 
essence of the ac Josephson effect 



Soon after these Josephson's predictions the micro- 
scopic theory of both dc [3] and ac [4,5] was constructed 
and these effects have been observed experimentally [7]. 
Huge number of publications as well as several mono- 
graphs are devoted to various aspects of these effects. It 
turned out that physics encoded in these phenomena is 
very rich and important for understanding of basic prop- 
erties of superconductivity itself. More than forty years 
after its discovery the Josephson effect still attracts at- 
tention of many researchers and keeps providing us with 
new interesting physics. 

In this paper I will discuss several new phenomena the- 
oretical understanding of which was achieved only very 
recently. In the next section I will very briefly review al- 
ready well known and established results which concern 
dc Josephson effect in various types of superconducting 
weak links. Sections II and III are devoted to possible 
new effects [8,9] which emerge and gain importance as 
one decreases the size of a weak link eventually turning 
it to a nanostructure with only few conducting channels. 
Fabrication of such quantum point contacts (QPC) - un- 
thinkable at the time of discovery of the Josephson effects 
- is now becoming a routine procedure. Hence, the new 
effects discussed here can be directly observed and inves- 
tigated in a modern experiment. 



I. INSTEAD OF INTRODUCTION 

Relatively soon after the Josephson's discovery it was 
understood that non-dissipative transport of Cooper 
pairs between two superconductors is possible not only 
through a (usually very thin) insulating barrier, but also 
in various other situations. One of such situations is re- 
alized in the so-called SNS structures, i.e. if a piece of a 
normal metal is placed in-between two superconductors. 
In contrast to tunnel junctions, in SNS systems at suf- 
ficiently low temperatures appreciable supercurrent can 
flow even though a normal layer can be as thick as few 
microns. This is because the wave function of Cooper 



1 



pairs or, more precisely, the anomalous Green function, 
penetrates into the normal metal from a superconduc- 
tor at the length ~ v F /T for ballistic and ~ y/D/T for 
diffusive metals (here and below D = vpl/3 and I are 
respectively the diffusion coefficient and the clastic mean 
free path). Clearly, at temperatures much lower than the 
critical temperature T c of a superconductor this length 
becomes large (as compared, e.g., to the superconducting 
coherence length), and macroscopic quantum coherence 
is established between two superconducting banks sepa- 
rated by a normal metal. 

Further studies revealed an interesting mechanism of 
Cooper pair transfer in such systems, ft turned out that 
the supercurrent flow is directly related to another fun- 
damentally important phenomenon: Andreev reflection 
[10]. Suffering Andreev reflections at both SN interfaces, 
quasiparticles with energies below the superconducting 
gap are effectively "trapped" inside the TV-layer and form 
a discrete set of levels [10]. It was demonstrated [11] that 
in the presence of the phase difference ip across the SNS 
junction these levels acquire a shift proportional to this 
phase difference. Thus, on one hand, the position of the 
quasiparticle energy levels in such systems can be tuned 
by passing the supercurrent and, on the other hand, the 
magnitude of this supercurrent can be established by tak- 
ing the derivative of the quasiparticle energy with respect 
to ip with subsequent summation over the whole energy 
spectrum. The microscopic theory [11,12] leads to the fol- 
lowing expression for the current density through clean 
SNS systems: 



J 



e 2 p 2 F v F 



if, —IT < Ifi < 7T. 



(3) 



67T 2 d 

This expression is valid at T — > and for TV-metal layers 
with thickness d S> £o ~ v F /A. The most important 
features of this result are (i) the strongly non-sinusoidal 
current-phase relation, cf. Eqs. (1) and (3) and (ii) the 
linear dependence of the current on the gap in the quasi- 
particle spectrum e qp <~ vp/d in the direction normal to 
NS interfaces. 

It is interesting that qualitatively both features (i) and 
(ii) survive not only for ballistic but also for diffusive 
SNS junctions even though in the latter case discrete 
Andreev levels are washed out due to elastic scattering 
of quasiparticles on impurities in the iV-metal. It was 
demonstrated microscopically [13-15] that at low tem- 
peratures T <C D/d 2 the current-phase relation in diffu- 
sive SNS junctions also deviates from the sinusoidal one 
[16] and the critical Josephson current is again propor- 
tional to the gap in the quasiparticle spectrum, in this 
case the Thouless energy e qp — D/d 2 . The exact value 
of the critical Josephson current in long diffusive SNS 
junctions can be established only numerically. One finds 
[15] 



10.82- 



eRN 



(4) 



where Rn is the junction normal state resistance. 

The above results - both for ballistic and diffusive lim- 
its - are valid for sufficiently long junctions. One can 
also decrease the thickness of the normal metal d and 
gradually crossover to the limit of short superconduct- 
ing constrictions. A microscopic description of the dc 
Josephson effect of such type of weak links was devel- 
oped by Kulik and Omel'yanchuk [17]. Also in such sys- 
tems at low temperatures the current-phase dependence 
deviates from sin ip and the critical current I C (T — ► 0) is 
again proportional to the combination e qp /eRN, where 
now e qp — A. A crossover between the two limits of long 
SNS junctions and short superconducting weak links can 
also be described microscopically. In the clean case this 
task can be trivially accomplished by solving the Eilcn- 
berger equations [18-20], while in the dirty limit one 
should make use of the Usadel equations [21] which can 
be solved only numerically. The latter task has recently 
been carried out in Ref. [15]. 

Let us also note that in all the above considerations 
inter-metallic interfaces were assumed to be perfectly 
transparent. It is also straightforward to generalize the 
analysis in order to account for electron scattering at the 
insulating barrier which can be present inside a weak 
link. For short superconducting junctions containing 
an insulating barrier with an arbitrary energy indepen- 
dent transmission the corresponding generalization has 
been worked out by Haberkorn et al. [22]. This anal- 
ysis yields a general formula for the Josephson current 
which matches with the Ambegaokar-Baratoff result [3] 
in the weak tunneling limit and crosses over to the Kulik- 
Omel'yanchuk's expression [17] for clean constrictions at 
transmissions approaching unity. It is interesting that 
the result [17] for diffusive constrictions can also be re- 
covered from the formula of Ref. [22] after its slight gen- 
eralization. In order to do so one should assume that 
the transmission is not the same for all conducting chan- 
nels but rather obeys the Dorokhov's distribution for- 
mula. Combining this formula with the expression [22] 
and summing over all conducting channels one arrives at 
the result [17] for diffusive weak links. 

One can also investigate transport properties of more 
complicated layered structures which contain both nor- 
mal metal layers and insulating barriers. For instance, 
SNS systems with one insulating barrier, such as SINS 
and SNINS were analyzed by a number of authors 
[23-29]. For an extended review summarizing various 
features of dc Josephson effect in different types of su- 
perconducting weak links and further references we refer 
the reader to Refs. [30-32]. 

Most of the results reviewed above were obtained al- 
ready long time ago and are by now well established 
and well understood. One can think that considering 
dc Josephson effect in even more complicated structures 
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like, for instance, SNS structures with two or three insu- 
lating barriers, may at most yield somewhat more cum- 
bersome expressions but would not allow to encover any 
new physics beyond what has already been understood 
in simpler situations. Below we will show that it is not 
so. Just on the contrary, in the next section we will 
demonstrate that qualitatively new effects may occur in 
SNS junctions with more than one insulating barriers, 
in particular provided the cross section of such junctions 
is reduced to be comparable to the square of the Fermi 
wavelength. 

II. JOSEPHSON EFFECT AND QUANTUM 
INTERFERENCE OF QUASIPARTICLES 

In this section we will analyze the dc Josephson effect 
in SNS systems which contain several insulating barri- 
ers. In this case electrons scattered at different barriers 
can interfere inside the junction. We will demonstrate 
that such interference may lead to qualitatively new ef- 
fects and cause severe modifications of the supercurrcnt 
across the junction. We will see that these modifications 
can go in both directions, i.e. the Josephson current can 
be dramatically decreased by destructive intereference of 
quasiparticles or, on the contrary, increased as a result 
of their constructive intereference. The first situation is 
realized for sufficiently short junctions, while for longer 
ones the second effect might become more pronounced. 

The phenomenon of quantum interference of quasipar- 
ticles is of primary importance for SNS structures with 
few conducting channels. The interest to such struc- 
tures grew considerably after several experimental groups 
[33-35] have succeeded in connecting a carbon nanotube 
to two superconductors and performing transport mea- 
surements in such systems. More conventional SNS 
structures with many conducting channels and several 
insulating barriers are also of considerable interest, for 
instance in relation to possible applications, see e.g. Ref. 
[36] and further references therein. We will demonstrate 
that for such systems quantum interference effects are 
also important provided there exist more than two scat- 
terers inside the junction. 

On a theoretical side a significant difficulty is that the 
powerful formalism of quasiclassical energy-integrated 
Eilcnberger Green functions [18-20,31] supplemented by 
the Zaitsev boundary conditions [37] cannot be directly 
applied to systems containing more than one insulating 
barrier. An important ingredient of the derivation [37] is 
the assumption that such barriers are located sufficiently 
far from each other, so that interference effects emerg- 
ing from electron scattering can be totally neglected. 
It is also essential that Zaitsev boundary conditions do 
not depend on the scattering phases. Since here we are 
just interested in investigating of quantum interference 
of quasiparticles we are not in a position to use the qua- 



siclassical Eilenberger formalism for our purposes. One 
possibility to circumvent this problem is to apply the 
formalism [38,39] within which the presence of an arbi- 
trary number of barriers in the system can be accounted 
for by linear boundary conditions. Another - even more 
straightforward - possibility to analyze the dc Josephson 
effect in structures with several insulating barriers is to 
directly solve the exact Gor'kov equations [40]. Here we 
will follow the second approach. 

The results presented in this section were obtained in 
collaboration with A.V. Galaktionov [8]. A similar ap- 
proach has also been used independently by Brinkman 
and Golubov [41]. 

A. General formalism 

In what follows we will assume that our system is uni- 
form along the directions parallel to the interfaces (coor- 
dinates y and z). Performing the Fourier transformation 
of the normal G and anomalous F + Green function with 
respect to these coordinates 

G Un {r,r') = J ^G w „(x,x',fc||)e lfc "^-i) 

we express the Gor'kov equations in the following stan- 
dard form 

/ iiO n — H A(x) \ / G Wn (x, x', fe|| ) \ = / 5(x - x') \ 

\ A*(x) iu> n + H c J \ F+ n (x,x', fe||) J \ J ■ 

(5) 

Here u n = (2n + l)nT is the Matsubara frequency, and 
A(x) is the superconducting order parameter. The Hamil- 
tonian H in Eq.(5) reads 

Here fc|| = fey — ep is Fermi energy, the term 

V(x) accounts for the external potentials (including the 
boundary potential). The Hamiltonian H c is obtained 
from H (6) by inverting the sign of the electron charge e. 

As usually, it is convenient to separate fast oscillations 
of the Green functions oc exp(±ik x x) from the envelope 
of these functions changing at much longer scales as com- 
pared to the atomic ones. Then one can construct a 
particular solution of the Gor'kov equations (5) in the 
following form 

^_ 2 {x)g 2 (x')e- lk * {x - x "> if x > x' (7) 

and 
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G Wn (a;,a;',fe||) 



Tp +2 (x)f 2 (x')e ik ^ x - x '^ iix<x'. (8) 



These functions satisfy Gor'kov equations at x ^ x' . 
Here Tp + arc two linearly independent solutions of the 
equation 



iuj n -Hl A(x) 
A*(x) iu n + H a ±c I r 



(9) 



The solution Tp +1 ) does not diverge at x — > +00, while 
Tp +2 is well-behaved at x — > —00. Similarly, two linearly 
independent solutions 2 do not diverge respectively 
at x — > —00 and x — ► +00. 
In eq. (9) we defined 



(10) 



Here k x — mv x — ^Jk 2 F — fejj, V^(x) represents a slowly 

varying part of the potential which does not include fast 
variations which may occur at metallic interfaces. The 
latter will be accounted for by the boundary conditions 
to be considered below. 

The functions f 1,2(2) and gi,2(x) are determined with 
the aid of the continuity condition for the Green functions 
at x = x' and the condition resulting from the integration 
of S(x — x') in eq.(5). 

A general solution of the Gor'kov equations has the 
form 



Gui n (x, x ) 
F+ n (x,x') 

[h(x')Tp +1 (x) 



Gu n (x, x ) 
F+ (x,x'\ , 

V ' / / part 

i(x')Tp +2 (x)]e 



■k x x 



+ 



(11) 



[h(x')<P-i(x) + k(x')f_ 2 (x)]e- 



ik x x 



For systems which consist of several metallic layers the 
particular solution is obtained with the aid of the pro- 
cedure outlined above provided both coordinates x and 
x' belong to the same layer. Should x and x' belong to 
different layers, the particular solution is zero because in 
that case the 5-function in eq. (5) fails. The functions 
h,2,i,i(x') in each layer should be derived from the proper 
boundary conditions. These are just the matching condi- 
tions for the wave functions on the left and on the right 
side of a potential barrier, respectively A\ exp(iki x x) + 
B\ exp(— iki x x) and A 2 enp(ik 2x x) + B 2 exp(— ik 2x x). is 
These conditions have the standard form (see e.g [42]): 



A 2 = aAt + pBu B 2 = (3*A X + a*B x , 

,|2 ^lx 



M 2 - 



k 



2x 



(12) 



The equations 



R 



D = l-R 



kla 



k2x\a\ 



(13) 



define respectively the reflection and transmission coeffi- 
cients of the barrier. Applying these boundary conditions 
at each insulating barrier one uniquely determines all the 
unknown functions in eq. (11) and thereby completes the 
construction the Green functions of our problem. For fur- 
ther details we refer the reader to Ref. [8]. 

We are now in a position to specify the general ex- 
pression for the Josephson current across ballistic SNS 
junctions which contain an arbitrary number of insulat- 
ing barriers. In what follows we will assume that a thin 
specularly reflecting insulating barriers (/) are situated 
at both SN interfaces. Additional such barriers can also 
be present inside the TV-metal. Transmissions of these 
barriers my take any value from zero to one. We also as- 
sume that electrons propagate ballistically between any 
two adjacent barriers and that no electron-electron or 
electron-phonon interactions are present in the normal 
metal. For simplicity we will restrict our attention to the 
case of identical superconducting electrodes with singlet 
isotropic pairing and neglect suppression of the supercon- 
ducting order parameter A in the electrodes close to the 
SN interface. The phase of the order parameter is set to 
be —tp/2 (+(p/2) in the left (right) electrode. As before, 
the thickness of the normal layer will be denoted by d. 

Employing the standard formula for the current den- 
sity 



le 
m 



d 2 k 



(2n)< 



(Vx' - Vx) x ^ x G Wn (x,x'M). 



(14) 



and making use of the expressions for the Green func- 
tions, one arrives at the following result 



J = 4eT^ 



w„>0 



k x dk x simp 
2ir cos ip + W' 



(15) 



where the function W depends on the number of insulat- 
ing barriers. This function will be specified below for the 
case of two and three barriers. 

Note that the integral over k x in eq. (15) can be re- 
placed by a sum over independent conducting channels 



2tt 



J 



N 



k x dk x (...) -> ^(...), 



(16) 



where A is the junction cross section. In this case Di i2 
and Ri,2 may also depend on the channel index m. 
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B. SINI'S junctions with few conducting channels 

Let us first consider SNS junctions with two insulat- 
ing barriers, one at each NS interface. In this case the 
function W in (15) takes the form 



W = - r4r cos x 



+ 



0^(1 + i2x)(l + R 2 ) 



D 1 D 2 A 2 
2(1 - R 1 R 2 ) n n iu r , 



+ 



DiD 2 



D X D 2 A 2 

lu^DiD 2 2u> n d 
— cosh 

v x 

. , 2uj n d 
sinn . 



(17) 



Here \ = 2k x d + (j) is the phase of the product a^fhctiPi- 
Eqs. (15), (17) provide a general expression for the dc 
Josephson current in SINI'S structures valid for arbi- 
trary transmissions D\ and D 2 . 

Let us first analyze the above result for the case of 
one conducting channel N — 1 . We observe that the first 
term in eq. (17) contains cos{2k x d+(j)) which oscillates at 
distances of the order of the Fermi wavelength. Provided 
at least one of the barriers is highly transparent and/or 
(for sufficiently long junctions d > £o) the temperature is 
high T 3> vp/d this oscillating term is unimportant and 
can be neglected. However, at lower transmissions of 
both barriers and for relatively short junctions d<vp/T 
this term turns out to be of the same order as the other 
contributions to W (17). In this case the supercurrent 
is sensitive to the exact positions of the discrete energy 
levels inside the junction which can in turn vary consid- 
erably if d changes at the atomic scales <~ 1/kp. Hence, 
one can expect sufficiently strong sample-to-sample fluc- 
tuations of the Josephson current even for junctions with 
nearly identical parameters. 

Let us first consider the limit of relatively short 
SINI'S junctions in which case we obtain 



e A T sin w 
I = — tanh 



V 



VA 
~2T 



(18) 



where we defined 



2%) 



1 - Tsin 2 (vV2) 



(19) 



and an effective normal transmission of the junction 

D X D 2 



T 



1 + R 1 R 2 + 2y/7TJl 2 ~cosx' 



(20) 



Eq. (18) has exactly the same functional form as the re- 
sult derived by Haberkorn et al. [22] for SIS junctions 
with an arbitrary transmission of the insulating barrier. 
This result is recovered from our eqs. (18), (20) if we as- 
sume e.g. Di <C D 2 in which case the total transmission 
(20) reduces to T ~ D 1 . 

As we have already discussed the total transmission T 
and, hence, the Josephson current fluctuate depending on 



the exact position of the bound states inside the junction. 
The resonant transmission is achieved for 2k x d+cf> = ±7r, 
in which case we get 



T 

1 X i- 



D X D 2 



(21) 



This equation demonstrates that for symmetric junctions 
D\ = D 2 at resonance the Josephson current does not 
depend on the barrier transmission at all. In this case 
T rcs = 1 and our result (18) coincides with the for- 
mula derived by Kulik and Omel'yanchuk [17] for bal- 
listic constrictions. In the limit of low transmissions 
Di,2 *C 1 we recover the standard Brcit-Wigner formula 
T ICS = \D\D 2 j(D\ + D 2 ) 2 and reproduce the result ob- 
tained by Glazman and Matveev [43] for the problem of 
resonant tunneling through a single Anderson impurity 
between two superconductors. 

Note that our results (18-20) also support the conclu- 
sion reached by Beenakker [44] that the Josephson cur- 
rent across sufficiently short junctions has a universal 
form and depends only on the total scattering matrix of 
the weak link which can be evaluated in the normal state. 
Although this conclusion is certainly correct in the limit 
d — > 0, its applicability range depends significantly on 
the physical nature of the scattering region. From eqs. 
(15), (17) we observe that the result (18), (19) applies 
at d <C Co not very close to the resonance. On the other 
hand, at resonance the above result is valid only under 
a more stringent condition d -C £o£>max, where we define 
Anax =max(£>i,L> 2 ). 

Now let us briefly analyze the opposite limit of suf- 
ficiently long junctions d » £o- Here we will restrict 
ourselves to the most interesting case T = 0. From eqs. 
(15), (17) we obtain 



I = 



ev x sin ip 
ndzi 



arctan y/z 2 /zi 



(22) 



\/z 2 jz\ 

21,2 - cos 2 (^/2) + -p-p- (R+ ± 2^R X R 2 cos(x)) , 

where R + = Ri + R 2 . For a fully transparent channel 
Di = D 2 = 1 the above expression reduces to the well 
known Ishii-Kulik result [12,11] 



ev x ip 
nd ' 



-7T < If < IT, 



(23) 



whereas if one transmission is small fli«l and D 2 w 1 
we reproduce the result [23] 



ev x D\ sin<£ 
2d ' 



(24) 



Provided the transmissions of both S'-interfaces are low 
Di.2 <C 1 we obtain in the off-resonant region 



-^ Dl D 2 sm V T[ X ], 



(25) 
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where T[x] is a 27r-periodic function defined as 

X 



T M = - 



sm x 



— 7T < % < 7T. 



(26) 



In the vicinity of the resonance ||x| — 7r j < D max the 
above result does not hold anymore. Exactly at reso- 
nance x = i 71 " we get 



ev x y/DiD 2 sinip 



1/2- 



(27) 



Ad { cos 2 f 



For a symmetric junction £>i j2 = -D this formula yields 



ev x D sin(<p/2) 
2d : 



-IT < tp < IT, 



(28) 



while in a strongly asymmetric case D\ <C D 2 we again 
arrive at the expression (24). This implies that at reso- 
nance the barrier with higher transmission D 2 becomes 
effectively transparent even if D 2 -C 1- We conclude that 
for Di t 2 <C 1 the maximum Josephson current is pro- 
portional to the product of transmissions D\D 2 off reso- 
nance, whereas exactly at resonance it is proportional to 
the lowest of two transmissions D\ or D 2 . 

We observe that both for short and long SINI'S junc- 
tions interference effects may enhance the Josephson ef- 
fect or partially suppress it depending on the exact po- 
sitions of the bound states inside the junction. We also 
note that in order to evaluate the supercurrent across 
SINI'S junctions it is in general not sufficient to de- 
rive the transmission probability for the corresponding 
NINI'N structure. Although the normal transmission 
of the above structure is given by eq. (20) for all values of 
d, the correct expression for the Josephson current can be 
recovered by combining eq. (20) with the results [22,44] 
in the limit of short junctions d <C D$o only. In this case 
one can neglect suppression of the anomalous Green func- 
tions inside the normal layer and, hence, the information 
about the normal transmission turns out to be sufficient. 
On the contrary, for longer junctions the decay of Cooper 
pair amplitudes inside the ./V-layer cannot be anymore 
disregarded. In this case the supercurrent will deviate 
from the form (18) even though the normal transmission 
of the junction (20) will remain unchanged. This devia- 
tion becomes particularly pronounced for long junctions, 
i.e. for d » ( out of resonance and for d 3> -D£o at 
resonance. 

Generalization of the above results to the case of an 
arbitrary number of independent conducting channels 
TV > 1 is trivial: The supercurrent is simply given by the 
sum of the contributions from all the channels. These 
contributions are in general not equal because the phase 
factors x = 2k x d + (j> change randomly for different chan- 
nels. Hence, mesoscopic fluctuations of the supercurrent 
should become smaller with increasing number of chan- 
nels and eventually disappear in the limit of large TV. 



In the latter limit the Josephson current is obtained 
by averaging over all values of the phase x- This limit 
was already studied in details [8,41] and will not be con- 
sidered here. We will only point out that - as it was 
demonstrated in Ref. [8] - in the limit N — > 00 interfer- 
ence effects are effectively averaged out and exactly the 
same result can be reproduced by means of the Eilcn- 
berger formalism supplemented by Zaitsev boundary con- 
ditions. We also worthwhile to emphasize that the latter 
statement applies only to the junctions with two insulat- 
ing barriers. Below we will show that for system with 
more than two barriers quasiparicle interference effects 
turn out to be even more significant, and the correct re- 
sult for the current cannot be recovered with the aid of 
Zaitsev boundary conditions even in the limit N — > 00. 



C. Josephson current in SINI' NI" S junctions 

Let us now turn to SNS structures with three insu- 
lating barriers. As before, two of them are located at 
SN interfaces, and the third barrier is inside the ./V-layer 
at a distance d\ and d 2 respectively from the left and 
right SN interfaces. The transmission and reflection co- 
efficients of this intermediate barrier are denoted as Do 
and i?o = f — Do, whereas the left and the right barri- 
ers are characterized respectively by D\ = 1 — R\ and 
D 2 = I-R2. 

The supercurrent is calculated along the same lines as 
it was done for the case of two barriers. The final result 
is again expressed by eq. (15), where the function W is 
now defined by a substantially more cumbersome expres- 
sion than one for the two barriers case. This expression 
was evaluated in Ref. [8] and will not be presented here. 
We will go over to the final results. 



1. One channel limit 

Let us first discuss the case of one conducting chan- 
nel. In the limit of short junctions d <C £oDmax we again 
reproduce the result (18) where the total effective trans- 
mission of the normal structure with three barriers takes 
the form 



2tit()t2 



1 + ht t 2 + C(<£>i,2,io,i,2) 



(29) 



where 



C = cosxi^l - tl){\ - t\) + cos X2 ^(l - tl)(l - t\) 

+ (cos xi cos X2 - t sin xi sin X2 ) \J {I - t\){l - t\). (30) 

Here we define i ,i,2 = A), 1,2/(1 + #0,1,2) and xi,2 = 
2k x di t 2 + <j>i,2- For later purposes let us also perform 



G 



averaging of this transmission over the phases xi,2- Wc 
obtain 



(V 



2tit()t2 



t 2 t 2 t 2 



(31) 



In particular, in the case of similar barriers with small 
transparencies -Do, 1,2 ~fCl the average normal trans- 
mission of our structure is (T) ~ D 3 ! 2 . Suppression of 
the average transmission below the value ~ D is a result 
of destructive interference and indicates the tendency of 
the system towards localization. 

Let us now proceed to the limit of a long junction 
di,2 ^> £0 an d T = 0. In the off-resonant region for 
d\ = d,2 we find 



ev x D 1 D D 2 sunp T[xi] - T[% 2 ] 



87rdi 



cos X2 - cosxi 



(32) 



This expression diverges at resonance (i.e. at Xi — n or 
X2 — ti") where it becomes inapplicable. In the resonant 
region % 2 ~ 7r we obtain 



ev x yJDiD D2 sin ip 



4^2(1 + cos Xl )(T-i-sin 2 (^/2)) 



(33) 



2. Many channel junctions 

As it was already discussed, in the many channel limit 
it is appropriate to average the current over the scat- 
tering phases. Practically in any realistic physical real- 
ization the widths d\ and d 2 fluctuate independently on 
the atomic scale. In this case averaging over \i and X2 
should also be performed independently If d\ and d 2 do 
not change on the atomic scale but are incommensurate, 
independent averaging over the two phases is to be per- 
formed as well. Independent averaging cannot be fulfilled 
only in (physically irrelevant) case of strictly commensu- 
rate d\ and rf 2 which will not be considered below. 

Technically independent averaging over the scattering 
phases X\ = x and \2 — Ax amounts to evaluating the 
integral of the expression l/[t + cosxcos(Ax)] from x = 
to some large value x = L. At A = 1 the result of this 
integration is L/^/t(l + 1). However, if A is irrational, 
the integral approaches the value 2LK(l/t 2 )/wt, where 
K{h) = F(n/2,h) is the complete elliptic integral. 

Let us assume that the transparencies of all three in- 
terfaces are small as compared to one. After averaging 
over the two scattering phases we arrive at the final ex- 
pression for the current 



eki „ . m A 2 T ^ 

w„>0 ™ 



rA 2 sin%/2) 



(34) 



D, 



eff 



/ ndny/D DiD 2 . (35) 
Jo 



Hence, for similar barriers we obtain the dependence 
J oc £> 3 / 2 rather than J oc D (as it would be the case 
for independent barriers). The latter dependence would 
follow from the calculation based on Zaitsev boundary 
conditions for the Eilcnberger propagators. We observe, 
therefore, that quantum interference effects decrease the 
Josephson current in systems with three insulating bar- 
riers. This is essentially quantum effect which cannot 
be recovered from Zaitsev boundary conditions even in 
the multichannel limit. This effect has exactly the same 
origin as a quantum suppression of the average normal 
transmission (T) due to localization effects. Further lim- 
iting expressions for short junctions can be directly re- 
covered from Eq. (31). 

We also note that the current-phase relation (34) devi- 
ates from a pure sinusoidal dependence even though all 
three transmissions are small -Do. 1,2 « 1 At T = the 
critical Josephson current is reached at ip ~ 1.7 which 
is slightly higher than ir/2. Although this deviation is 
quantitatively not very significant, it is nevertheless im- 
portant as yet one more indication of quantum interfer- 
ence of electrons inside the junction. 

Finally, let us turn to the limit of long junctions rfi ;2 » 
Co- We again restrict ourselves to the case of low trans- 
parent interfaces. At high temperatures T 3> vp/2ndi i 2 

d 

we get J oc D a DiD 2 e «< T > , where d — d\ + d 2 and 
£(T) = vf/{2ttT). In this case the anomalous Green 
function strongly decays deep in the normal layer. Hence, 
interference effects are not important and the interfaces 
can be considered as independent from each other. In the 
opposite limit T <C Dvp/d, however, interference effects 
become important, and the current becomes proportional 
to D 5 / 2 rather than to D 3 . Explicitly, at T -> with the 
logarithmic accuracy we get 



J 



ekpVp sin ip 

\§TT 2 \Jd\d,2 



[ dfi^D^y^lnDQ 1 . (36) 
^0 



where wc define the effective transmission 



We see that, in contrast to short junctions, in the limit 
of thick normal layers interference effects increase the 
Josephson current as compared to the case of indepen- 
dent barriers. The result (36), as well as one of Eqs. (34) 
(35) cannot be obtained from the Eilenberger approach 
supplemented by Zaitsev boundary conditions. 



D. Some conclusions 

By directly solving the Gor'kov equations we evaluated 
the dc Josephson current in SNS junctions containing 
two and three insulating barriers with arbitrary trans- 
missions, respectively SINI'S and SINI'NT'S junc- 
tions. Our results can be directly applied both to the 
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junctions with few conducting channels (such as, e.g., 
superconductor-carbon nanotube-superconductor junc- 
tions [33-35]) and to more conventional SNS structures 
in the many channel limit. We have demonstrated that 
an interplay between the proximity effect and quantum 
interference of quasiparticlcs may play a crucial role in 
such systems causing strong modifications of the Joseph- 
son current. 

For the system with two barriers and few conducting 
channels we found strong fluctuations of the Josephson 
critical current depending on the exact position of the 
resonant level inside the junction. For short junctions 
d <C £oD at resonance the Josephson current does not 
depend on the barrier transmission D and is given by 
the standard Kulik-Omel'yanchuk formula [17] derived 
for ballistic weak links. In the limit of long SNS junc- 
tions d 3> Co resonant effects may also lead to strong 
enhancement of the supercurrent, in this case at T — > 
and at resonance the Josephson current is proportional 
to D and not to D 2 as it would be in the absence of 
interference effects. 

While the above results for few conducting channels 
cannot be obtained by means of the approach employing 
Zaitsev boundary conditions, in the many channel limit 
and for junctions with two barriers the latter approach 
does allow to recover correct results. This is because the 
contributions sensitive to the scattering phase are effec- 
tively averaged out during summation over conducting 
channels. 

Quantum interference effects turn out to be even more 
important in the proximity systems which contain three 
insulating barriers. In this case the quasiclassical ap- 
proach based on Zaitsev boundary conditions fails even 
in the limit of many conducting channels. In that limit 
the Josephson current is decreased for short junctions 
(J cx D 3 / 2 ) as compared to the case of independent bar- 
riers (J oc D). This effect is caused by destructive inter- 
ference of electrons reflected from different barriers and 
indicates the tendency of the system towards localization. 
In contrast, for long SNS junctions with three barriers 
an interplay between quantum interference and proxim- 
ity effect leads to enhancement of the Josephson current 
at T — > 0: We obtained the dependence J cx D 5 / 2 instead 
of J cx D 3 for independent barriers. 

III. PARITY AFFECTED JOSEPHSON 
CURRENT 

Let us now turn to a different issue which - to the 
best of our knowledge - was not yet attracted attention 
in the literature. Namely, we will discuss an interplay 
between the parity effect and the dc Josephson current 
in superconducting weak links. The results presented in 
this section have been obtained in collaboration with S.V. 
Sharov [9]. 



It is well known that thermodynamic properties of iso- 
lated superconducting systems are sensitive to the parity 
of the total number of electrons [45,46] even though this 
number J\f is macroscopically large. This parity effect is 
a direct consequence of the fundamental property of a su- 
perconducting ground state described by the condensate 
of Cooper pairs. The number of electrons forming this 
condensate is necessarily even, hence, for odd Af at least 
one electron always remains unpaired having an extra 
energy equal to the superconducting energy gap A. At 
sufficiently low temperatures a clear difference between 
the superconducting states with even and odd Af was 
demonstrated experimentally [46,47]. 

Can the supercurrent be affected by this parity effect? 
At the first sight the answer to this question should be 
negative because of the fundamental uncertainty relation 
5N8ip > 1. Should the electron number Af be fixed, fluc- 
tuations of the superconducting phase <p become large 
disrupting the supercurrent in the system. On the other 
hand, suppressing fluctuations of the phase ip will destroy 
the parity effect because of large fluctuations of Af. 

Despite that, below we will demonstrate that in certain 
superconducting structures the parity effect can coexist 
with the non- vanishing supercurrent. Consider a super- 
conducting system which can support circular persistent 
currents (PC). For An example is provided by an iso- 
lated superconducting ring pierced by the magnetic flux 
<f> in which case circulating PC is induced in the ring. 
In accordance with the number-phase uncertainty rela- 
tion the global superconducting phase of the ring fluc- 
tuates strongly in this case, however these fluctuations 
are decoupled from the supercurrent and therefore can 
be integrated out without any influence on the latter. In 
what follows we will show that the parity effect may sub- 
stantially modify PC in superconducting nanorings, in 
particular for odd number of electrons. 

A. Parity projection formalism 

In order to systematically investigate the influence of 
the electron parity number on persistent currents in su- 
perconducting nanorings we will employ the well known 
parity projection formalism [49-51]. Recapitulating the 
key points of this approach we will closely follow Ref. 
[50]. 

The grand canonical partition function Z(T,fi) = 
Tre -/3(«-AiA0 ig ii n k e d to the canonical one Z(T,Af) by 
means of the following equation 

Z(T lA1 ) = f>(T,A0cx P (^Y (37) 

Here and below TL is the system Hamiltonian, Af is the to- 
tal number of electrons and (3 = 1/T. Inverting this rela- 
tion and defining the canonical partition functions Z e and 
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Z Q respectively for even (TV = N e ) and odd {N = Af ) 
ensembles, one gets 

7T 

Z e /o(T) = ^ J due-^>° u Z e/o (T,iTu), (38) 



where the upper/lower sign corresponds to the even/odd 
ensemble and we have defined 



'e/o 



0fi, 



/o 



A»(*x) 



/an 



f/b 



f*(*x) 



where 



Z e/o (T, M ) = l'&{[l±(-l)^|, 



-(3(H-nN) 



B. Parity effect in nanorings and blocking of the 
supercurrent 



-(Z(T,/i)±Z(T,/i + *7rT)) 



(39) 



are the parity projected grand canonical partition func- 
tions. For J\f 1 it is sufHcient to evaluate the integral in 
(38) within the saddle point approximation which yields 



Z e /o(T) ~ e -P$i*io-»eioX*io) ^ 



(40) 



where £l e / = —ThiZ e / (T, /x) arc the parity projected 
thermodynamic potentials. They can be presented in the 
form 



O e/o = 0/-Tln 



±.(i± e -m*-n f )^ 



where 



n //6 = -TIn 



Tr 



{(±1) 



M p -0{H-tiM) 



(41) 



(42) 



Chemical potentials /i e / Q are defined by the saddle point 
condition N e/o = -d9. e / {T,^ e/o )/d^ e / . 

The main advantage of the above analysis is that it 
allows to express the canonical partition functions and 
thermodynamical potentials in terms of the parity pro- 
jected grand canonical ones thereby enormously simpli- 
fying the whole calculation. We further note that f2/ is 
just the standard grand canonical thermodynamic poten- 
tial and Qb represents the corresponding potential linked 
to the partition function Z(T, \i + iirT). It is easy to see 
[50] that in order to recover this function one can evalu- 
ate the true grand canonical partition function Z(T,/j,), 
express the result as a sum over the Fermi Matsubara 
frequencies ojj = 2-KT(m + 1/2) and then substitute the 
Bose Matsubara frequencies uib = 2irTm instead of the 
Fermi ones. This procedure will automatically yield the 
correct expression for Z(T,fi + iirT) and, hence, for f2 b . 

Having found the thermodynamic potentials for the 
even and odd ensembles one can easily determine the 
equilibrium current /. Here we will be interested in de- 
scribing the currents flowing in isolated superconducting 
rings pierced by the external magnetic flux <fr x . Then 
in the case of even/odd total number of electrons one 
obtains 



h/o - !f ± 



I 



f 



e P(si b -n f ) ± i 



(43) 



Let us now make use of the above general expressions 
and investigate the influence of the parity effect on PC 
in superconducting nanorings with quantum point con- 
tacts (QPC). Before turning to concrete calculations we 
shall specify the model for our system. We shall con- 
sider mesoscopic superconducting rings with cross sec- 
tion s and perimeter L — 2irR. The rings will be as- 
sumed sufficiently thin, i.e. \fs <C \l, where Al is the 
London penetration length. Superconductivity will be 
described within the (parity projected) mean field BCS 
theory. At sufficiently low temperatures this description 
is justified provided quantum phase slips (QPS) [52-54] 
in nanorings can be neglected. This requirement in turn 
implies that the ring cross section should be sufficiently 
large. With the aid of the results [52] one concludes 
that the QPS tunneling amplitude remains exponentially 
small provided the condition s 3> X% \J £o Jl is satisfied. 
Here Xf is the Fermi wavelength, £ ~ «f/A is the co- 
herence length and I is the electron elastic mean free path 
assumed to be shorter than £o- For generic systems QPS 
effects can usually be neglected provided the transversal 
size of the wire/ring y/~s exceeds ~ 10 nm. Hence, the to- 
tal number of conducting channels in the ring 7V r ~ s/\p 
should inevitably be large N T ^$> 1. In addition, the ring 
perimeter L should not be too large, so that one could 
disregard the QPS-induced reduction of the PC ampli- 
tude [54]. Finally, we will neglect the difference between 
the mean field values of the BCS order parameter for 
the even and odd ensembles [49,50]. This is legitimate 
provided the volume of a superconducting ring is large 
enough, V = Ls 3> 1/^A, where v is the density of states 
at the Fermi level and A is the BCS order parameter for 
a bulk superconductor at T = 0. All these requirements 
can easily be met in a modern experiment. 

The task at hand is now to evaluate the thermody- 
namic potentials Qf/b- Within the mean field treatment 
these quantities can be expressed in terms of the excita- 
tion energies and the superconducting order parameter 
A(r). One finds [50] 



flf = fi-2T^]ln(2cosh— ) , 

k 

n b = n-2T]Tln (2sinh||) , 



(44) 
(45) 
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where O = / d 3 r|A(r)| 2 /.g+Tr{£}, g [ s the BCS coupling 
constant and £ is the single-particle energy operator: 



1 

2^i 



■ih 



d_ 

dr 



A(r)\ +U(r)- l i, 



(46) 



A(r) is the vector potential and U(r) describes the po- 
tential profile due to disorder and interfaces. 
The excitation spectrum Sk has the form 



£fc = e(p) 



+ + a 2 , 



(47) 



where p is a quasiparticle momentum, £ = (p 2 — ft) /2m, 
and /i = /ti($ x ) — mvs 2 /2. The superconducting velocity 
vector «s is oriented in the direction along the ring and 
is defined by the well known expression 



v s = 



2mR 



5x 

$0 



(48) 



This expression as well as the excitation spectrum (47) 
are the periodic functions of the flux $> x with the period 
equal to the superconducting flux quantum $o = hc/2e. 

Consider the most interesting case T — > 0. Making use 
of the above expressions one easily finds 

I e = ev s g e s, n e = H(Q e ) + mv 2 s /2 (49) 

for the even ensemble and 

Vf 1 I Vs I 
Io = ev s g s - e— sgn(v s ), Q a = Q + rr • (50) 

Li V Dp 

for the odd one. Here g e / a = N e / /V are the electron 
densities for the even and odd ensembles, g is the grand 
canonical electron density at T = 0, = ^J2[ijm and 
vs is assumed to be small as compared to the critical 
velocity vc = A /pp. We also note that the second Eq. 
(50) is an implicit equation for the chemical potential /i . 

Eq. (49) - being combined with (48) - coincides with 
that obtained for the grand canonical ensemble. In par- 
ticular, the current I e represents the well known "saw 
tooth" dependence on magnetic flux. In contrast, for 
odd ensembles there exists an additional flux-dependent 
contribution to PC (50) which cannot be viewed just as 
a renormalization of g a . 

Unfortunately this parity effect is rather small in mul- 
tichannel rings [16]. Estimating the leading contribution 
to I e / a as I <~ evpNr/L, wc find 

(I e ~ I )/I ~ l/N T < 1. 

The results (49)-(50) hold as long asT« hvp/L. At 
higher temperatures the parity effect gets even smaller 
and eventually disappears at temperatures exceeding the 
parameter [46] T* w A/ \n(vV\/AT*). The correspond- 
ing expressions are readily obtained within our formal- 
ism, but we will not consider them here. 



Rather we turn to a somewhat different system - a su- 
perconducting ring interrupted by QPC - in which the 
parity effect turns out to play a much more important 
role. In this case the thermodynamic potential of the 
system fi consists of two different contributions [55] 



(51) 



respectively from the bulk part of the ring and from QPC. 
The optimal value of the phase difference (p across QPC 
is fixed by the condition dSl/dip — which reads 



—c- 



d<P x 



h dip 



(52) 



Here we made use of the fact that the thermodynamic 
potentials of the ring depend both on & x and ip only via 
the superfluid velocity vs = (l/47rmi?)((/?— 27r<j> :E /<I>o)> in 
which case one can put d/d$ x = —(2e/Tic)(d/dtp). The 
left-hand side of Eq. (52) represents the current flowing 
inside the superconducting ring 1^ = —cdfl^d^ x ~ 
(evpN T /L)((p — 27r$ :E /$o)- This value should be equal to 
the current across QPC which is given by the right-hand 
side of Eq. (52). Estimating the maximum value of the 
latter for a single channel QPC as 2eTA/h, we obtain 



p~2ir-^, if L«L 
<£>o 



If ~ 27m, 



if L > L* , 



(53) 
(54) 



where L* = ^oN T /T ^> £o- In a more general case of 
QPC with N conducting channels in the expression for 
L* one should set T — ► J2n ^n- 

In what follows we will consider the most interesting 
limit N < iV r and L < L* . Due to Eq. (53) in this 
case the dependence I e /o(®x) is fully determined by the 
current-phase relation for QPS which can be found by 
means of Eq. (43) with Ifj b = — (2e/h)dfl^J h /dip. It is 
convenient to employ the formula [8] 



N 



l f/b 



sin</3 



2e 



& i ^ cos 99 + W n (uj f / b ) 

n=l Uf/b 



(55) 



In the case of short QPS one has W n (oj) = (2/T n )(l + 
w 2 /A 2 ) — l, where T n is the transmission of the n-th con- 
ducting channel. Substituting this function into (55) and 
summing over uij one recovers the standard result [17,22] 



N 



dip 



2T 



where 



s n (<p) = AJ 1 - T n sin 2 (<p/2). 



(56) 



(57) 



The same summation over Bose Matsubara frequencies 
ujb yields 
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N 



h = 



2e d£ n&) cQth £n(y) 



h 



n=l 



2T 



(58) 



Finally, the difference £7^ — 0/ = is evaluated as 
a sum of the ring and QPS contributions. 

The latter is found by integrating If/^p) over the phase 
difference p: 



N 



= 2T^]lncoth 



71=1 



2T 



(59) 



while the former is defined by the standard expression 
[50] 



7 2tt 



P lnfcoth £(p) 



Combining all these results with Eq. (43) we get 



2e ^— \ ut, 



N 



d£n ^ tanh- £ " ( ^ 



1± 



(coth : 



dp 



>(v)\2 



2T 



2T 



e mi » r f n(coth 

i=l 



2T 



) 2 ±1 



(60) 



Eq. (60) represents the central result of this section. 
Together with Eq. (53) it establishes the complete de- 
pendence of PC on the magnetic flux $ in isolated su- 
perconducting nanorings with QPC. 

Consider the most interesting limit T — > 0. In this case 
for the even number of electrons in the ring PC is given 
by the expression (56) which coincides identically with 
that for grand canonical ensembles [17,22]. On the other 
hand, for the odd number of electrons PC will acquire 
an additional contribution which turns out to be most 
important for the case of single channel QPS N = 1. In 
that case the expression in the square brackets of Eq. 
(60) reduces to zero, i.e. PC will be totally blocked by 
the odd electron. Thus, we predict a novel mesoscopic 
effect - parity affected blocking of PC in superconducting 
nanorings with QPC. 

This result has a transparent physical interpretation. 
Indeed, it is well known [48] that the result (56) can be 
expressed via the contributions of discrete Andreev levels 
E±(p) = ±AV(p) inside QPS as 



2e 



dp 



f-(E- 



(61) 



where T>(p) is defined in Eq. (19). Using the Fermi filling 
factors for these levels f±{E±) = [1 + exp(£; ± (^)/T)]- 1 
one arrives at Eq. (56). If we now fix the number of 
electrons inside the ring and consider the limit T — > 
the filling factors will be modified as follows. For the 



even number AT all electrons are paired occupying states 
with energies below the Fermi level. In this case one has 
/_(£?_) = 1, f + (E + ) = 0, the current is entirely de- 
termined by the contribution of the quasiparticle state 
75 _ and Eq. (61) yields the same result as one for the 
grand canonical ensemble. By contrast, in the case of 
odd number of electrons one electron always remains un- 
paired and occupies the lowest available energy state - 
in our case E + - above the Fermi level. Hence, for odd 
Af one has f±(E±) = 1, the contributions of two quasi- 
particle energy states in Eq. (61) exactly cancel each 
other, and the current across QPS remains zero for any 
p or the magnetic flux & x . This is just the blocking ef- 
fect which we have already obtained above from a more 
formal consideration. 

For TV > 1 and/or at non-zero temperatures this 
parity-affected blocking of PC becomes incomplete. But 
also in this case the parity effect remains essential at 
temperatures T < T* substantially affecting, e.g., the 
current-phase relation for QPC. For T > this relation 
will deviate from the grand canonical one both for even 
and odd ensembles [9]. 

Finally, we turn to superconducting rings containing a 
piece of a normal metal. Here we only restrict our at- 
tention to transparent SNS junctions with length of the 
normal metal d 3> £o- I n this case for uo <C A one has 
Wi(u) — cosh(2ujd/vF)- Substituting this function into 
(55) and repeating the the whole calculation as above, in 
the limit T^Owe obtain 



evpN 
nd 



evpN 
nd 



7TSgn</3\ 

N ) 



(62) 



These results apply for — 7r < ip < ir and should be 2n- 
periodically continued otherwise. We observe that the 
current I e again coincides with that for the grand canon- 
ical ensembles [11], while for odd Af the current-phase re- 
lation is shifted by the value ir/N. This shift has a simple 
interpretation as being related to the odd electron contri- 
bution (2e/h)dE /dp from the lowest (above the Fermi 
level) Andreev state E (p) inside the SNS junction. Un- 
like in QPC, this contribution does not compensate for 
the current from other quasiparticle states. Rather it pro- 
vides a possibility for a parity-induced 7r-junction state 
in our system. Indeed, according to Eq. (62) for single 
mode SNS junctions the "saw tooth" current-phase rela- 
tion will be shifted exactly by ir. For more than one con- 
ducting channel N > 1 within the interval — ir < p < tt 
there exists a twofold degenerate minimum energy (zero 
current) state [56] occurring at p = ±ir/N. In the par- 
ticular case N = 2 the current-phase relation I {p) turns 
7r-periodic. 

The well known feature of superconducting rings in- 
terrupted by a 7r-junction is the possibility to develop 
spontaneous supercurrent in the ground state [57]. Al- 
though this feature is inherent to any type of 7r-junctions, 
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in the case of the standard sinusoidal current-phase rela- 
tion such spontaneous supercurrents can occur only for 
sufficiently large values of the ring inductance C [57]. In 
contrast, in the situation studied here the spontaneous 
current state is realized for any inductance of the ring 
because of the non-sinusoidal dependence I (<fi) (62). 

In order to demonstrate that let us assume that no ex- 
ternal flux is applied to our system. Then at T —* the 
energy of an SNS ring with odd number of electrons can 
be written in the form 



E n = 



$ 2 irhv F N 



2cC 



$ sgn<I> 
2N 



(63) 



where $ is the flux related to the circular current flowing 
in the ring. Minimizing this energy with respect to $, 
one easily observes that a non-zero spontaneous current 



± 



evF 



1 + 



2ev F N C 



d $ 



(64) 



should flow in the ground state of our system. This is yet 
one more remarkable consequence of the parity effect: 
Just by changing TV from even to odd one can induce 
non-zero PC without any external flux <fr x . In the limit 
of small inductances C <C $od/evpN - which is easy to 
reach in the systems under consideration - the value of 
/ does not depend on the number of channels N and is 
given by the universal formula / = ±evp/d. For generic 
parameters this value can easily be as large as I ~ 10 na. 

In summary, new physical effects emerge from an inter- 
play between the electron parity number and persistent 
currents in superconducting nanorings. These effects can 
be directly tested in modern experiments and possibly 
used for engineering new types of superconducting flux- 
charge qubits. 
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